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Abstract 

Supersymmetric extensions of the Standard Model when combined with the Peccei-Quinn 
solution to the strong CP problem necessarily contain also the axino, the fermionic partner 
of the axion. In contrast to the neutralino and the gravitino, the axino mass is generically 
not of the order of the supersymmetry breaking scale and can be much smaller. The axino is 
therefore an intriguing candidate for a stable superpartner. The axinos are a natural candidate 
for cold dark matter in the Universe when they are generated non-thermally through out-of- 
equilibrium neutralino decays or via a competing thermal production mechanism through 
scatterings and decays of particles in the plasma. We identify axino masses in the range of 
tens of MeV to several GeV (depending on the scenario) as corresponding to cold axino relics 
if the reheating temperature Tr is less than about 5 x 10 4 GeV. At higher Tr and lower mass, 
axinos could constitute warm dark matter. In the scenario with axinos as stable relics the 
gravitino problem finds a natural solution. The lightest superpartner of the Standard Model 
spectrum will remain stable in high-energy detectors but may be either neutral or charged. 
The usual constraint Q x /i 2 < 1 on the relic abundance of the lightest neutralino does not hold. 



CERN-TH/2001-059 
February 2001 

1 Based on invited plenary talks at the 4th International Workshop on Particle Physics and the Early 
Universe (COSMO-2000), Cheju, Cheju Island, Korea, September 4-8, 2000 and the 3rd International 
Workshop on the Identification of Dark Matter (IDM-2000), York, England, 18-22 September 2000. 



1 Introduction 



The nature of dark matter (DM) in the Universe remains unknown. Its relic abundance 
is probably in the range 0.1 <Q x h 2 < 0.3 or so and is most likely "cold". It presumably 
consists for the most part of some weakly interacting massive particles (WIMPs). 

WIMPs do not necessarily have to interact via weak interactions per se. Very much 
weaker interaction strengths can also lead to interesting cosmological abundances. One 
expects that they should preferably be electrically and color neutral, and therefore be 
non-baryonic, otherwise they would dissipate their kinetic energy. 

Supersymmetry (SUSY) provides a natural context for WIMPs. In the presence 
of i?-parity, the (massive) lightest supersymmetric particle (LSP) is stable and may 
substantially contribute to the relic mass density. It may therefore constitute the DM in 
the Universe. Among SUSY WIMPs, the neutralino remains the most popular choice by 
being perhaps the most "genuine" WIMP. The gravitino is another interesting candidate 
although it generically suffers from the (in)famous "gravitino problem". 

Recently, it has been pointed out|], that another well- motivated SUSY particle 
can be a plausible cold DM (CDM) candidate. The particle in question is the axino, 
a fermionic SUSY partner of the axion. Axionic particles arise naturally in various 
superstring models and, historically first, in the Peccei-Quinn (PQ) solution to the 
strong CP problem. 

It is also well-justified to consider the axino as the LSP since its mass is basically 
a free parameter which can only be determined in specific models. Since current LEP 
bounds on the neutralino are in excess of 32 GeV, or so, it becomes increasingly plausible 
that there may well be another SUSY particle which is lighter than the neutralino, and 
therefore a candidate for the LSP and dark matter. The axino is just such a candidate. 

In this talk I will first describe the main properties of the new candidate. Next I will 
discuss ways of producing axinos in the early Universe and the expected abundance. 
Finally, I will make some remarks about the ensuing implications for cosmology and 
phenomenology. 

2 Axino 

Axino is a superpartner of the axion pseudoscalar. It is thus a neutral, i?-parity negative, 
Majorana chiral fermion. There exist several SUSY and supergravity implementations 
of the well-known original axion models (KSVZ[[3| and DFSZ[Q). (Axion/axino-type 
supermultiplets also arise in superstring models.) In studying cosmological properties of 
axinos, we will concentrate on KSVZ and DFSZ-type models where the global U(l) PQ 
symmetry is spontaneously broken at the PQ scale f a . A combination of astrophysical 
(white dwarfs, etc) and cosmological bounds gives 10 9 GeV < f a ~ 10 12 GeV[|] although 
the upper bound can be significantly relaxed if inflation followed the decoupling of 
primordial axionic particles and the reheating temperature Tr <C f a . 
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The two main parameters of interest to us are the axino mass and coupling. The 
mass rria strongly depends on an underlying model and can span a wide range, from 
very small (~ eV) to large (~ GeV) values. What is worth stressing is that, in contrast 
to the neutralino and the gravitino, axino mass does not have to be of the order of the 
SUSY breaking scale in the visible sector, M SU sy ~ 100 GeV — 1 TeV|B|, [7], |8j. The basic 
argument goes as follows. In the case of unbroken SUSY, all members of the axion 
supermultiplet remain degenerate and equal to the tiny mass of the axion given by the 
QCD anomaly. Once SUSY is softly broken, superpartners acquire mass terms. Since 
the axino is a chiral fermion, one cannot write a dimension-four soft mass term for the 
axino. (For the same reason there are no soft terms for, e.g., the MSSM higgsinos.) The 
lowest-order term one can write will be a non-renormalizable term of dimension-five. 
The axino mass will then be of order M| USY // a ~ lkeV[§, f|. 

However, in specific SUSY models there are normally additional sources of axino 
mass which can generate much larger contributions to m„ at one-loop or even the tree- 
level. One-loop terms will always contribute but will typically be < Msusy (KSVZ) 
or even <C Msusy (DFSZ). Furthermore, in non- minimal models where the axino mass 
eigenstate comes from more than one superfield, m„ arises even at the tree-level. In 
this case m s can be of order Msusy but can also be much smaller. In a study of 
cosmological properties of axinos, it therefore makes sense to treat their mass as a 
basically free parameter. 

Axino couplings to other particles are generically suppressed by l/f a - For our pur- 
pose the most important coupling will be that of axino-gaugino-gauge boson interactions 
which can written as a dimension-five term in the Lagrangian 



where B denotes the bino, the fermionic partner of the U(l)y gauge boson B, which is 
one of the components of the neutralino, g stands for the gluino and N — 1(6) for the 
KSVZ (DFSZ) model. One can show that the SU{2)l coupling can be rotated away 
so long as one discusses cosmological properties of the axino at large temperatures. 
Depending on a model, one can also think of terms involving dimension-four operators 
coming, e.g., from the effective superpotential where ^> is one of MSSM matter 

(super)fields. However, axino production processes coming from such terms will be 
suppressed at high energies relative to processes involving Eq. (Jl|) by a factor m%/s 
where s is the square of the center of mass energy. We will comment on the role of 
dimension-four operators again below but, for the most part, mostly concentrate on the 
processes involving axino interactions with gauginos and gauge bosons, Eq. ([]]), which 
are both model- independent and dominant. 
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3 Axino Production 



Particles like axinos and gravitinos are somewhat special in the sense that their interac- 
tions with other particles are very strongly suppressed compared to the SM interaction 
strengths. Therefore such particles remain in thermal equilibrium only at very high 
temperatures. In the particular case of axinos (and axions), their initial thermal popu- 
lations decouple at [[/J] Td ~ 10 10 GeV. At such high temperatures, the axino co-moving 
number density is the same as the one of photons. In other words, such primordial 
axinos freeze out as relativistic particles. Rajagopal, Turner and Wilczek (RTW)[[/j] 
pointed out that, in the absence of a subsequent period of inflation, the requirement 
that the axino energy density is not too large (Qa ^ 1) leads to m„ < 2 keV and the cor- 
responding axinos would be light and would provide warm or even hot dark matter 0. 
We will not consider this case in the following primarily because we are interested in 
cold DM axinos. We will therefore assume that the initial population of axinos (and 
other relics, like gravitinos) present in the early Universe was subsequently diluted away 
by an intervening inflationary stage and that the reheating temperature after inflation 
was smaller than Ti> It also had to be less than f a , otherwise the PQ would have been 
restored thus leading to the well-known domain wall problem associated with global 
symmetries. 

In order to generate large enough abundance of axinos, one needs to repopulate the 
Universe with them. There are two generic ways of doing achieving this. First, they can 
be generated through thermal production, namely via two-body scattering and decay 
processes of ordinary particles and sparticles still in thermal bath. (Despite the name 
of the process, the resulting axinos will typically be already out-of-equilibrium because 
of their exceedingly tiny couplings to ordinary matter.) Second, axinos may also be 
produced in decay processes of particles which themselves are out-of-equilibrium. Such 
particles could for example be ordinary superpartners, the gravitino or the inflaton field. 
(By "ordinary" we mean a particle or its superpartner carrying only Standard Model 
quantum numbers.) Below we will concentrate on the first possibility in the case the 
decaying neutralino. 

3.1 Thermal Production 

After inflation the Universe can be re-populated with axinos (and gravitinos) through 
scattering and decay processes involving superpartners in thermal bath. As long as the 
axino co-moving number density ria is much smaller than n 7 , the number density of 
photons in thermal equilibrium, its time evolution will be adequately described by the 
Boltzmann equation 

—jj- + 3Hn d = + J « H Hci)^j + X^( r (* )) n i- ( 2 ) 



3 



Here H is the Hubble parameter, H(T) = J T , where g* is the effective relativistic 

degrees of freedom, cr(i + j —* 5 + - • •) is the scattering cross section for particles i, j into 
final states involving axinos, v IC \ is their relative velocity, Hi is the ith particle number 
density in thermal bath, Y(i — > a + • • •) is the decay width of the ith particle and (• ■ ■) 
stands for thermal averaging. (Averaging over initial spins and summing over final spins 
is understood.) Note that on the r.h.s. we have neglected inverse processes since they 
are suppressed by 

The main axino production channels are the scatterings of (s)particles described by 
a dimension- five axino-gaugino-gauge boson term in the Lagrangian ([!]). Because of the 
relative strength of a s , the most important contributions will come from 2-body strongly 
interacting processes into final states involving axinos, % + j —* a + • • •. (Scattering 
processes involving electroweak interactions are suppressed by both the strength of the 
coupling and a smaller number of production channels @.) The cross section can be 
written as 

a 3 

a n (s) = oCr n (s) (3) 

where y'i is the center of mass energy and n = A, . . . ,K refers to different channels 
which are listed in Table I in Ref . || . The diagrams listed in the Table are analogous to 
those involving gravitino production and we use the same classification. This analogy 
should not be surprising since both particles are neutral Majorana superpartners. 

In addition to scattering processes, axinos can also be produced through decays of 
heavier superpartners in thermal plasma. At temperatures T > m~ these are dominated 
by the decays of gluinos into LSP axinos and gluons. The relevant decay width is given 
by 

, a 2 ™£ ( m ?\ 3 

T(g a ^a + g b )= 5 ab ^^ a —* 1 - (4) 

1 } 1287T 3 (f a /N) y m|J K) 

and one should sum over the color index a, b = 1, • • • , 8. At lower temperatures m x < 
Tr < m~, neutralino decays to axinos also contribute while at higher temperatures they 
are sub-dominant. They only become important when the axino yield 



TP 

rr = (5) 

where s = (27r 2 /45)g s *T 3 is the entropy density, and normally g s * = g* in the early 
Universe, becomes too small to be cosmologically interesting. 

The results are presented in Fig. 1 for representative values of f a = 10 11 GeV and 
m~ = m~ = 1 TeV. The respective contributions due to scattering as well as gluino and 
neutralino decays are marked by dashed, dash-dotted and dotted lines. 

It is clear that at high enough Tr, much above m~ and m~, scattering processes 
involving such particles dominate the axino production. For Tr ^> nig, nig, y scat grows 
linearly as Tr becomes larger. In contrast, the decay contribution above the gluino mass 
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Figure 1: Y? p as a function of Tr for representative values of f a = 10 11 GeV and 
m~ = m~ = 1 TeV. 

threshold, Y dec ~ 5xl0~ 4 (Mpr^/m?), remains independent of Tr. At Tr roughly below 
the mass of the squarks and gluinos, their thermal population starts to become strongly 
suppressed by the Boltzmann factor e~ mj/T , hence a distinct knee in the scattering 
contribution in Fig. 1. It is in this region that gluino decays (dash-dotted line) given by 
Eq. ([|) become dominant before they also become suppressed by the Boltzmann factor 
due to the gluino mass. For m x < T R < rrig, nig, the axino yield is well approximated by 
F TP w Y dec ~5x 10~ 4 (MpT g /T^) e~ m s/ TR , and depends sensitively on the reheating 
temperature. 

At still lower temperatures the population of strongly interacting sparticles becomes 
so tiny that at Tr ~ m x neutralino decays start playing some role before they too become 
Boltzmann factor suppressed. We indicate this by plotting in Fig. 1 the contribution of 
the lightest neutralino (dotted line). It is clear that the values of Y" a TP in this region are 
so small that, as we will see later, they will not play any role in further discussion. We 
therefore do not present the effect of the decay of the heavier neutralinos. Furthermore, 
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model-dependent dimension-four operators will change axino production cross section 
at lower T R ~ M SU sy but will be suppressed at high temperatures. We have not studied 
this point yet. 

We emphasize that axinos produced in this way are already out of equilibrium. Their 
number density is very much smaller than n 7 (except Tr ~ 10 9 GeV and above) and 
cross sections for axino re-annihilation into other particles are greatly suppressed. This 
is why in Eq. (^) we have neglected such processes. Nevertheless, even though axinos 
never reach equilibrium, their number density may be large enough to give f2„ ~ 1 for 
large enough axino masses ( keV to GeV range) as we will see later. 



3.2 Non-Thermal Production 

The mechanism for non-thermal production (NTP) that we will consider works as fol- 
lows. Consider some lightest ordinary superpartner (LOSP). Because axino LSP cou- 
plings to everything else are suppressed by 1/ f a , as the Universe cools down, all heavier 
SUSY partners will first cascade-decay to the LOSP. The LOSPs then freeze out of 
thermal equilibrium and subsequently decay into axinos. 

A natural (although not unique), candidate for the LOSP is the lightest neutralino. 
For example, in models employing full unification of superpartner masses (like the 
CMSSM/mSUGRA), a mechanism of radiative electroweak symmetry breaking typi- 
cally implies fi 2 3> M 2 , where M\ is the bino mass parameter. As a result, the bino-like 
neutralino often emerges as the lightest ordinary superpartner || |l(], 11 . 



In the following we will assume that LOSP is the neutralino. It can decay to the 
axino and photon x — > 07 with the rate|], 3 

a 2 C 2 m 3 / m?\ 3 

r( x -> 57) = a ^° a f m * _ i-^f . (6) 

Here a em is the electromagnetic coupling strength, C axi = (C a yy/ cos 9w)Zu, with 
Zu standing for the bino part of the lightest neutralino. (We use the basis x% = 
ZiiB + Z i2 W 3 + Z i3 H® + Z i4L Ht (i = 1,2,3,4) of the respective fermionic partners 
(denoted by a tilde) of the electrically neutral gauge bosons B and W3, and the MSSM 
Higgs bosons and H t .) 

The corresponding lifetime can be written as 

,,2 \ -2 / r /m \ 2 /mnn„v\ 3 / ™2\- 3 



0.33 sec / a 2 em \ ( f a /N \ f 100 GeV 



HIT. 



^^> = MU7^J iraj {—) l 1 "^ ' < 7 > 

For large enough neutralino mass, an additional decay channel into axino and Z 
boson opens up but is always subdominant relative to x ~ > because of both the 
phase-space suppression and the additional factor of tan 2 ^^. As a result, even at 
m x 3> mz,rria, r(% — > aZ) ~ 3.35 r(x — > 07). It is also clear that the neutralino 
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lifetime rapidly decreases with its mass (~ 1/m 3 ). On the other hand, if the neutralino 
is not mostly a bino, its decay will be suppressed by the Z n - factor in C axi . 

Other decay channels are the decay into axino and Standard Model fermion pairs 
through virtual photon or Z but they are negligible compared with the previous ones. 
We will discuss them later since, for a low neutralino mass, i.e., long lifetime, they can, 
even if subdominant, produce dangerous hadronic showers during and after nucleosyn- 
thesis. 

Additionally, in the DFSZ type of models, there exists an additional Higgs-higgsino- 
axino couplings, which may open up other channels ||. These are model-dependent and 
I will not discuss them here. 

3.3 Constraints 

Several nontrivial conditions have to be satisfied in order for axinos to be a viable CDM 
candidate. First, we will expect their relic abundance to be large enough, Vtah 2 ~ 0.2. 
This obvious condition will have strong impact on other bounds. Next, the axinos 
generated through both TP and NTP will in most cases be initially relativistic. We 
will therefore require that they become non-relativistic, or cold, much before the era of 
matter dominance. Furthermore, since NTP axinos will be produced near the time of 
BBN, we will require that they do not contribute too much relativistic energy density to 
radiation during BBN. Finally, axino production associated decay products will often 
result in electromagnetic and hadronic showers which, if too large, would cause too 
much destruction of light elements. In deriving all of these conditions, except for the 
first one, the lifetime of the parent LOSP will be of crucial importance. 

A detailed discussion of the bounds would take too much time, and space, that is 
available. I will therefore merely summarize the relevant results. First, the condition 
that the axinos give a dominant contribution to the matter density at the present time 
can be expressed as mgYa ~ 0.72 eV (f2 5 /i 2 /0.2) which applies to both TP and NTP 
relics. It is worth mentioning here that, for the initial population of axinos, the yield at 
decoupling is approximately Y & ~ ~ 2x 10 -3 which gives m s ~ 0.36 keV (flah 2 /0.2). 
This is an updated value for the RTW bound. Next, we want to determine the tem- 
perature of the Universe at which the axinos will become non-relativistic. In nearly all 
cases axinos are initially relativistic and, due to expansion, become non-relativistic at 
some later epoch which depends on their mass and production mechanism. In the case 
of TP, axinos are not in thermal equilibrium but, since they are produced in kinetic 
equilibrium with the thermal bath, their momenta will have a thermal spectrum. They 
will become non-relativistic when the thermal bath temperature reaches the axino mass, 

NTP axinos generated through out-of-equilibrium neutralino decays will be produced 
basically monochromatically, all with the same energy roughly given by m x /2, unless 
they are nearly mass-degenerate with the neutralinos. This is so because the neutralinos, 
when they decay, are themselves already non-relativistic. Thus, due to momentum red- 
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shift, axinos will become non-relativistic only at a later time, when p(Tnr) — mj. The 
temperature T NR can be expressed as|Q| 

T NR = 4.2 x lO^maCaYvZu (m x /100GeV) 1/2 (lO 11 GeV// a /iV 

This epoch has to be compared to the matter-radiation equality epoch given by 
T cq = 1.1 eV (Qah 2 /0.2) which holds for both thermal and non-thermal production. In 
the TP case one can easily see that T NR > T eq is satisfied for any interesting range of 
mj. In the case of NTP the condition Tnr T cq is satisfied for 

-»™ ^ra' /2 (i^)(¥)- oo 

If axinos were lighter than the bound (^), then the point of radiation- matter equality 
would be shifted to a later time around Tnr. Note that in this case axino would not 
constitute cold, but warm or hot dark matter. In the NTP case discussed here other 
constraints would however require the axino mass to be larger than the above bound, 
so that we can discard this possibility. 

BBN predictions provide further important constraints on axinos as relics. In the 
case of non-thermal production most axinos will be produced only shortly before nucle- 
osynthesis and, being still relativistic, may dump too much to the energy density during 
the formation of light elements. In order not to affect the Universe's expansion during 
BBN, axino contribution to the energy density should satisfy p&/ Pv < 5N V , where p v 
is the energy density of one neutrino species. Agreement with observations of light 
elements requires 8N V = 0.2 — 1. This leads top 

-> isivv 1 1 f lOOGeV V'V fJN \ ( Sl^\ 
m 5 181keV^^-(^— j ^^—j. (10) 

Finally, photons and quark-pairs produced in NTP decays of neutralinos, if produced 
during or after BBN, may lead to a significant depletion of primordial elements. One of- 
ten applies a crude constraint that the lifetime should be less than about 1 second which 
in our case would provide a lower bound on m x . First, photons produced in reaction 
X — > cff carry a large amount of energy, roughly m x /2. If the decay takes place before 
BBN, the photon will rapidly thermalize via multiple scatterings from background elec- 
trons and positrons. The process will be particularly efficient at plasma temperatures 
above 1 MeV which is the threshold for background ee pair annihilation, and which, 
incidentally, coincides with time of about 1 second. But a closer examination! 12j shows 



that also scattering with the high energy tail of the CMBR thermalize photons very 
efficiently and so the decay lifetime into photons can be as large as 10 4 sec. By com- 
paring this with Eq. (|7|) we find that, in the gaugino regime, this can be easily satisfied 
for m x < mz- It is only in a nearly pure higgsino case and mass of tens of GeV that 
the bound would become constraining. We are not interested in such light higgsinos for 
other reasons, as will be explained later. 
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A much more stringent constraint comes from considering hadronic showers from qq- 
pairs. These will be produced through a virtual photon and Z exchange, and, above the 
kinematic threshold for x —* a>Z, also through the exchange of a real Z-boson. Here the 
discussion is somewhat more involved and the resulting constraint strongly depends on 
m x . One can show0 that at the end one gets roughly m a > 360 MeV for m x < 60 GeV 
which gives the strongest bound so far. However, the bound on m„ decreases nearly 
linearly with m x and disappears completely for m x > 150 GeV. 

In summary, a lower bound m„ > (9(300 keV) arises from either requiring the axinos 
to be cold at the time of matter dominance or that they do not contribute too much to 
the relativistic energy density during BBN. The constraint from hadronic destruction 
of light elements can be as strong as m s > 360 MeV (in the relatively light bino case) 
but it is highly model-dependent and disappears for larger m x . 

3.4 Relic Abundance from Thermal and Non-Thermal Pro- 
duction 

In the TP case the axino yield is primarily determined by the reheating temperature. 
For large enough Tr (Tr ^> m~, m~), it is proportional to T^/f 2 . In contrast, the NTP 
axino yield is for the most part independent of Tr (so long as Tr ^> If, the neutralino 
freezeout temperature). In the NTP case, the yield of axinos is just the same as that of 
the decaying neutralinos. This leads to|]] Vtah 2 = ma/m x Q x h 2 , where Q x h 2 stands for 
the abundance that the neutralinos would have had today had they not decayed into 
axinos. 

In order to be able to compare both production mechanisms, we will therefore fix the 
neutralino mass at some typical value. Furthermore we will map out a cosmologically 
interesting range of axino masses for which fi^ TP ~ 1. Our results are presented in 
Fig. 2 in the case of a nearly pure bino. We also fix m x = 100 GeV and f a = 10 11 GeV. 
The dark region is derived in the following way. It is well known that Q x h 2 can take 
a wide range of values spanning several orders of magnitude. In the framework of the 
MSSM, which we have adopted, global scans give Vt x h 2 < 10 4 in the bino region at 
m x < 100 GeV. (This limit decreases roughly linearly (on a log- log scale) down to 
~ 10 3 at m x ~ 400 GeV.) For m x = 100 GeV we find that the expectation tt^ TP h 2 ~ 1 
gives 

10 MeV <rria< m x . (11) 

We note, however, that the upper bound Q x h 2 < 10 4 comes from allowing very large 
Msusy (i.e., sfermion and heavy Higgs masses) in the range of tens of TeV. Restricting 
all SUSY mass parameter below about 1 TeV reduces Q x h 2 below 10 2 and, accordingly, 
increases the lower bound m„ > 1 GeV. Still, for the sake of generality, we have the 
much more generous bound ([11]) in Fig. 2 but marked a low range of m„ with a light 
grey band to indicate the above point. 
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Likewise, for reheating temperatures just above Tf standard estimates of Q x h 2 be- 
come questionable. We have therefore indicated this range of Tr with again a light grey 
color. It has also been recently pointed out ||13|| that even in the case of very low reheat- 
ing temperatures Tr below the LOSP freezeout temperature, a significant population of 
them will be generated during the reheating phase. Such LOSPs would then also decay 
into axinos as above. We have not considered such cases in our analysis and accordingly 
left the region Tr < Tf blank even though in principle we would expect some sizeable 
range of f2 a /i 2 there. 

We can see that for large Tr, the TP mechanism is more important than the NTP 
one, as expected. Note also that in the TP case the cosmologically favored region 
(0.2 <Qah 2 ^0.4) would form a very narrow strip (not indicated in Fig. 2) just below 
the Qj p = 1-boundary. In contrast, the NTP mechanism can give the cosmologically 
interesting range of the axino relic abundance for a relatively wide range of m a so long 
as Tr < 5 x 10 4 GeV. Perhaps in this sense, the NTP mechanism can be considered as 
somewhat more robust. 

3.5 Conclusions 

The intriguing possibility that the axino is the LSP and dark matter possesses a number 
of very distinct features which makes this case very different from those of both the 
neutralino and the gravitino. In particular, the axino can be a cold DM WIMP for a 
rather wide range of masses in the MeV to GeV range and for relatively low reheating 
temperatures Tr < 5 x 10 4 GeV. As Tr increases, thermal production of axinos starts 
dominating over non-thermal production and the axino typically becomes a warm DM 
relic with a mass broadly in a keV range. In contrast, the neutralino is typically a cold 
DM WIMP. 

Low reheating temperatures would favor baryogenesis at the electroweak scale. It 
would also alleviate the nagging "gravitino problem". If additionally it is the axino that 
is the LSP and the gravitino is the NLSP, the gravitino problem is resolved altogether 
for both low and high Tr. 

Phenomenologically, one faces a well-justified possibility that the bound Q x h 2 < 1, 
which is often imposed in constraining a SUSY parameter space, may be readily avoided. 
In fact, the range Q x h 2 ^> 1 (and with it typically large masses of superpartners) 
would now be favored if axino is to be a dominant component of DM in the Universe. 
Furthermore, the lightest ordinary superpartner could either be neutral or charged but 
would appear stable in collider searches. 

The axino, with its exceedingly tiny coupling to other matter, will be a real challenge 
to experimentalist. It is much more plausible that a supersymmetric particle and the 
axion will be found first. Unless the neutralino (or some other WIMP) is detected in 
DM searches, the axino will remain an attractive and robust candidate for solving the 
outstanding puzzle of the nature of dark matter in the Universe. 
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Figure 2: The thick solid line gives the upper bound from thermal production on the 
reheating temperature as a function of the axino mass. The dark region is the region 
where non-thermal production can give cosmologically interesting results (Q^ TP h 2 ~ 1) 
as explained in the text. We assume a bino-like neutralino with m x = 100 GeV and 
f a = 10 11 GeV. The region of T R > T f is somewhat uncertain and has been denoted 
with light-grey color. A sizeable abundance of neutralinos (and therefore axinos) is 
expected also for T R < T f but has not been calculated. The vertical light-grey band 
indicates that a low range of m a corresponds to allowing SM superpartner masses in the 
multi-TeV range, as discussed in the text. Division of hot, warm and cold dark matter 
by the axino mass shown in lower left part is for axinos from non-thermal production. 
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